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THEORY  OF  THE  Z-MODE  PROPAGATION 


Notation 


The  following  frequencies  will  appear  in  our  formulae, 


f  «  operating  wave  frequency  (in  Me)  = 


■■N 


Ne 

nra 


1/2 


3  1/2  3 

9  X  10  (N)  »  plasma  frequency  (in  Me  if  N  «  el/cm  ) 


eH 


^  =  2,8  Hq  ■  cyclotron  frequency  (in  Me  if  H^  in  Gauss) 


H  2»tmc 
V  B  collision  frequency 


=  2  7tfj^  '  *  critical  collision  frequency. 


We  will  consider  a  plane  wave  propagating  in  the  plane  of  the  magnetic 
meridian  in  a  horizontally  stratified  ionosphere. 

The  angle  of  incidence  (Qg) ,  the  angle  of  the  wave  normal  in  the  layer  of 
index  of  refraction  n  (6) ,  and  the  angle  of  the  magnetic  field  (a)  in  the  selected 
co-ordinate  system  are  shown  in  Fig,  1. 


Figure  1 


The  following  compact  symbols  will  also  be  used. 


X  -  .  _Nei 

f2  5Cmf2 


Y  ■  Jl  a  -I-  9.  ,  Y,  ■  Ycosoi  ,  Yx  ■  Ysintx 
f  2jlmcf 


U  «  1-iZ  •  1-i  V  ^  (i.e.,  Z  -  y  J 

IWl  ’  zjtf 

COS0JJ  *  C,  sinSg  *  S 


Development  of  the  basic  equations 


The  equation  of  motion  of  a  plane  wave  in  a  magnetoionic  medium  is 


Introducing  the  volume  polarization  P 


Ne?^  we  obtain 


r* 


imVcu 

“TIT" 


P  a  eE  + 


e  icu 
c"  Me 


and,  therefore 


or 


-moj^  ,,  j  ^  N  Tt  .  imU)^  1 
^  45tNe2  ^  47tNe2  «*>  me 


-X 

Tji 


T  •  \ip  -  iY^? 


writing  eqn.  [2]  in  matrix  form  we  have 


X 

47r 


u 

-i^L 

iY-p 


iY^  -iYT 


U 

0 


0 

u 


7  2 

and  inverting  the  matrix  equation,  provided  ^  Y  ,  we  obtain: 


-X 

4  7X  U(u2-y2) 


-iY^U  iYjU 
iY^U  u2.y2  -YlYt 

-iYjU  -YlYj  U^-yJ 


The  matrix  eqn,  [3]  can  be  written  as: 


-4 


where  matrix  a  is  called  the  susceptibility  matrix.  We  also  have; 

D  =  6.E  =  E  +  40tP  *  E  ♦  4^t.o<E  ■  (1  ♦  43tot)  E  and  thus; 

€  s  1  +  4ato<  [4] 

where  €  is  now  the  complex  dielectric  tensor. 

From  eqns,  [3]  and  [4]  we  obtain  the  expression  for  £  ; 


- 

1  * 

iXY^U 

XYt-U 

^xx 

^xy 

^XZ 

U(U2-Y2) 

UCU'iY^) 

U(UiY2) 

£ 

c 

.  XYlU 

,  X(u2.y|) 

XYjYl 

'^yx 

eyy 

tyz 

U(u2-Y2) 

U{u2.y2) 

U(u2-y2) 

iXY^U 

XYtYl 

X(U^-yJ) 

^zx 

czy 

^zz 

u(u^-y2) 

u(u2-y2) 

^  ’  U(u2-y2) 

The  dependence  of  the  wave  functions  for  propagation  in  the  y-z  plane  through 
a  horizontally  stratified  ionosphere  (along  the  z  axis)  is; 

F  -  Fq  e  , 

where 

,  271  27tf  w 

h  «  -r-  ■  — r-  •  -r 
>  c  c 

so  that  we  have; 

^  “  i<*> »  37  "  0,  ^  -  -ikS. 
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Let  us  also  introduce  the  operator  -  =  q 

IK  dz 


Maxwell's  equations  in  the  absence  of  source  currents  and  changes  become: 


VaE  ■  -  i  li  .  .  H  [6] 

VaH  -  i  D  -  liiL  D  .  iidcE  [7] 

0  c  c 

Using  the  operators  defined  above,  eqns.  [6]  and  [7]  take  the  form  of  eqns,  [8]  and 
[9]  respectivel/o 

By  A  SE^  . 

-q  E,  .  -  Hy  [81 


-q  Hy  ♦  SHjj  »  -  ^xx^x  ■  ^xy^y  ”  ^xz^z 
-q  Hj5  »  ^yx^x  ^yy^y  *  ^yz^z 
SHx  •  ^zx^x  *  ^zyEy  + 

Eliminating  and  from  eqns.  [8]  and  [9]  we  obtain  a  system  of  four 


equations: 
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Substituting  from  eqn,  [5]  the  £•  tensor  elements  in  eqn,  [10]  and  writing 
it  also  in  a  matrix  form  we  have 


For  the  case  of  vertical  incidence  i.e,,  S  »  0,  C  «  1,  eqn,  [11]  takes  the  form: 
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— 

— 

“  — 

q 

0 

0 

1 

-E 

X 

0 

q 

1 

0 

iXYj^(U-X) 

U (U-X) ^-Y^ (U-X) 

n 

u 

U(u2-y2)-X(u2-Yl) 

10  2  2 
U(U^-Y-^)-X(U  -Y^) 

q 

u 

n 

X 

U(U-X)^-Uy2  +  XY^ 

-iXYL(U-X) 

0 

H 

V 

U(u2-y2).X(U^-Yl) 

U(U^-Y^)-X(U^-Y^) 

q 

[12] 


Finally  we  can  eliminate  H„  and  H  from  eqn. 


[12]  and  we  have: 


2  U(U-X)^-UY^  ♦  XY^ 

^  U(U^-Y2)-X(U^-Yf) 

XYl (U-X) 

U(U^-y2)-X(u2-Yl) 


XY^CU-X) 

^  2  9  2  2 

U(U  -Y‘^)-X(U^-Yp 

2  U(U-X)2-y2(U-X) 

U(U^-Y^)-X(U^-Y,^) 

L 


[13] 


Eqnsc  [11]  and  [13]  are  respectively  the  basic  equations  governing  the  propagation 
of  electromagnetic  waves  in  magnetoionic  media  for  oblique  and  vertical  incidence. 


The  Booker’s  quartic  and  its  solutions 


For  a  homogeneous  region  where  X,  Y,  U  are  constants,  the  z  dependence  is  only 
in  the  exponent  of  the  wave  function,  i.e,. 


8- 


-ikncosez 

‘IP  '  TIT  Hz  pQ® 


•ikncos6 

-ik 


>ikncos6z 


-ikncos6 


■  ncosdo  Fq® 


ncosQF 


and  therefore 


q  ■  ncos6. 

Requiring  the  determinant  of  eqn,  [11]  to  be  zero  we  obtain  a  quartic  equation  in  q 
called  the  Booker's  quartic. 


F(q) 


Ofq'^  +  ^q^  ♦  ^q^  +  Sq  +€ 


[14] 


where 

O  =  U(U^-Y^)-X(U^-Y^)  «  U^(U-X)-Uy2  +  XY? 

L  t. 

^  «  2SXYlY^ 

5  »  -2U[(cVx)  (U-X)-C^Y^]-X[y2-S^y2  + 

S  «  -2SC^XYj^Y^  » 

€  -  (C^U-X)  [(C^U-X)(U-X)-C^y2]-S2c2y|x. 

For  vertical  incidence  (S  ■  0  and  C  *  1) ,  ^and  6  become  zero  and  thus  we  have: 

F(q)  *0<q'^+^q^  +  €.  »0  [15] 


<X  .  U^(U-X)-UY^  +  XY^ 

5  »  -2U[(U-X)^-Y^]-X(Y^  +  Y^) 


where 
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e  .  (u-x) [(u-x)^-y2] 

Eqnc  [15]  is  also  obtained  by  equating  to  zero  the  determinant  of  eqn.  [13],  Thus, 

4  2 

for  vertical  incidence,  Booker's  quartic  becomes  a  biquadratic  ofq  *  ^  ^  ^  "0 

which  is  easily  solved  to  give  the  Appleton-Hartree  equation; 


1  - 


U  - 


2(U-X)  “  y4(U-X)2 


♦  y: 


[16] 


where  now,  since  cos6  «  1,  q^  »  n^  *  the  complex  index  of  refraction. 

When  collisions  are  neglected  (V  ) ,  so  that  U  »  1,  the  solutions  of  the 
Appleton-Hartree  equation  for  Y  ^  1,  are; 


n  «  0  at  X  ■  1  (ordinary) 
n«0  at  X»l-Y  (extraordinary) 


we  also  have  an  infinity,  n 


CD 


at  X  * 


_I_,  for  the  extraordinary, 

v2 


These 


solutions  hold  for  all  values  of  the  angle  o<  ,  with  one  exception,  the  case  when 
CX  =  0,  i,ec ,  Y^  •  0 


n 


2 


1 


X 

TTT 


and  the  roots  are; 

r*0  at  X-l+Y  (ordinary) 

1.  »0  at  X»l-Y  (extraordinary) 
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and  infinity  (n  «  oo)  at  X  *  1  for  the  extraordinary. 

We  see, therefore,  that  at  vertical  incidence,  the  extraordinary  wave  will 

always  be  reflected  (for  Y  ^  1)  at  X  ■  1  -  Y.  An  exponentially  decaying  part  of 

1  -  Y^ 

it  (evanescent  wave)  can  advance  farther,  but  it  will  be  absorbed  at  X  ■  . . . 

where  n  oo, 

The  ordinary  wave  in  the  case  of  vertical  incidence  will  always  be  reflected 
at  X  =  1  except  when  the  magnetic  field  is  also  vertical,  in  which  case  it  will  be 
reflected  at  X  ■  1  +  Y,  For  the  more  general  case  of  oblique  incidence  in  the  plane 
of  the  magnetic  meridian  the  four  roots  of  Booker's  quartic  correspond  to  the  four 
waves: 

1)  Ascending  ordinary 

2)  Descending  ordinary 

3)  Ascending  extraordinary 

4)  Descending  extraordinary 

When  two  of  the  roots  are  equal,  we  have  either  reflection,  e.g.,  roots 
of  ascending  and  descending  ordinary  waves  equal,  or  coupling  e.g.,  roots  of 
ascending  ordinary  and  extraordinary  equal. 

To  have  two  roots  of  eqn,  [14]  equal  we  must  also  have 


In  order  to  satisfy  eqns,  [14]  and  [17]  simultaneously  we  must  have 
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8a^  -3^^  60<S-P^ 


16«e  -p5 


6aS-p^  4a€,-^^  +  2p5 


16(xe-pS 


8€^  -3^^ 


The  solutions  of  this  determinant  have  been  investigated  by  M.  L.  V,  Pitteway  (1959), 
In  the  simpler  case  of  vertical  incidence  we  simply  have  ps.S  =  o,  and  the  above 


determinant  reduces  to: 


Cxe(g^-40(6.)^  «  0, 


which  is  satisfied  when: 


«  s  0:  (l-X)-Y^  ♦  XY?  -  0  i,e.,  X  -  i-lJL  (Infinity) 

“ —  TT77 

6=0:  (1-X) [(l-X)^-Y^]  »  0  i.e,,  X«l,  X«l-Y,  X-l+Y  (Reflections) 

i~^ - 

2  -4a6  «  0:  X(U-X)  \ ,  — ...T . ♦  Y?  ■  0  which  at  X  ■  1  yields: 

-  V4(U-X)2  b 


^*27“ 


(Coupling) 


Diagrammatic  representation 


We  will  now  investigate  the  variation  of  q  with  X  in  eqns,  [14]  and  [15] 


3«-4  .  ^3  .  3X  2  d^  dfe 


=  0  at  the  point  where  F(q)  has  a  double  root.  At  that  point,  ^  becomes  infinite 


-12- 


and  the  q-X  curves  have  a  vertical  tangent.  If,  however,  the  numerator  of  eqn,  [18] 

is  zero  also  at  the  same  point,  ^  is  not  infinite  and  the  tangent  is  not  vertical 

dX 

F(q),  neglecting  collisions,  can  be  factored  at  X  =  1  to  give; 


y 

F(q)  «  (q^-c2)(q-S  -  0. 

Coupling  will  take  place  at  the  double  root,  i.e,,  at  q  ■  S  Inserting  in  the 

Yj 

numerator. of  eqn,  [18]  the  above  value  of  q  and  the  derivatives; 


2SYlYt. 


rdocl 

x«i 

L35rJx«i 

[^]x=i  '  *  s2y2.c2y2 

fdSl 
bxJx-1 

Pdel 
LhtJx.i 


-2SC"’Yj^Y^ 


S'^  +  C^Y^  -  S^C^Y^, 


We  obtain; 


sin^<X 


which  becomes  equal  to  zero  when; 


S  *  sin0^  ■ 
o 


*  \f  ^  si 
-  Vy  +  1 


sino^ . 


Thus,  at  X  *  1,  for  Y  ^  1  and  sin0Q  s  ^  sinCX ,  we  have  complete  coupling 

instead  of  reflection.  No  collisions  are  necessary  for  this  coupling*  For  the 


[19] 
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case  of  vertical  incidence  (6^  ■  0)  complete  coupling  without  collisions  occurs  only 
when  0<  =  0,  When  oi  ^  0,  coupling  occurs  at  X  =  1  only  in  the  presence  of  the 
required  collision  frequency 


V 


2  3ff, 


H 


sin^ 

2cos(X 


No  such  simple  formula  exists  for  the  case  of  oblique  incidence  when  sinS^ 

In  Fig.  2  we  will  plot  the  curves  of  q  vs,  X  for  vertical  and  oblique  incidence, 
illustrating  schematically  the  regions  of  reflection  and  coupling.  The  curves  are 
derived  from  Booker's  quartic,  neglecting  the  collision  frequency,  i.e.,  U  ■  1, 

For  the  case  of  vertical  incidence  (6o  ■  0)  in  Fig,  2,  we  show  the  following 
situations? 


(a)  The  magnetic  field  makes  a  rather  large  angle  («)  with  the  vertical. 
Reflections?  Ordinary  at  X  =  1,  Extraordinary  at  X  »  1  -  Y, 

(b)  Angle  ”®  ^®'^®  coupling  between  the  ordinary  and 

extraordinary  waves,  producing  the  Z-mode  which  is  reflected  at 
X  »  1  +  Y,  The  shaded  areas  show  the  regions  of  coupling, 

(c)  Angle  (X  ■  0,  We  no  longer  have  the  ordinary  wave  reflected  at 
X  »  1, 


Note  that  the  vertical  incidence  curves  are  symmetric  as  the  ascending  and  descending 
waves  travel  the  same  path.  For  the  case  of  oblique  incidence  we  have 
correspondingly  ? 

(d)  TTie  dashed  lines  represent  the  case  where  Bq  is  much  larger  than  «  , 

The  reflections  are  no  longer  atX»l-Y,  1,  and  1  ♦  Y  but  at 


si 
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lower  values  of  X,  given  respectively  to  a  first  approximation  by 
X  «  (1  -  Y)cos6q,  COS0O,  and  (1  +  Y)coseoo  The  solid  lines 
represent  the  case  where  sinSjj  ^  J  ^  sino<  .  Here  6q  is  quite 
small  and  the  curves  approach  those  of  vertical  incidence »  For  all 
the  values  of  0  sin"^|^  ^  ^  ^  sincxj  ,  the  ordinary  is 

reflected  at  X  ■  1« 

(e)  Here  sine^  ^  ^  sino<  and  we  have  coupling,  ioe,,  part  of  the 

ordinary  wave  is  reflected  near  X  *  1  and  part  of  the  wave  couples 

to  the  extraordinary  which  propagates  farther  and  is  finally  reflected 
near  X  *  (1  +  Y)coseg, 

(f)  This  is  the  case  where  sinS^  ■  \j^  ^  ^  sinoco  At  X  «  1  we  now  have 
complete  coupling  instead  of  reflection.  We  notice  that  the 
oblique  incidence  curves  are  not  symmetric  because  the  descending 
path  is  different  from  the  ascending  one. 

When  a  moderate  number  of  collisions  are  present,  —  0.01,  the  curves 

are  very  similar  to  those  shown  with  the  exception  of  the  coupling  regions  where  a 
full  wave  treatment  is  necessary. 

The  coupling  parameters  have  been  determined  for  the  case  of  vertical  incidence 
through  the  full  wave  treatment  of  the  problem,  but  no  similar  analysis  has  yet  been 
accomplished  for  the  case  of  oblique  incidence  where  the  physical  mechanism  is 
analogous  but  the  mathematics  are  very  much  more  complex. 
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Coupling 

We  will  now  examine  the  coupoing  process  for  the  case  of  vertical  incidence. 
An  ordinary  wave  incident  normally  on  an  interface  separating  two  regions  of 
different  indices  of  refraction  (n_  »  1  -  ,  ,X  1  and  accordingly  different 
polarizations  (p^) »  will  have  to  produce  an  extraordinary  wave  in  order  to  satisfy 
the  boundary  conditions, 

The  extraordinary  wave  produced  will  increase  as  the  change  of  ordinary 
polarization  (Apo)  increases,  or  as  the  difference  between  the  ordinary  (p^)  and 
extraordinary  (p^)  polarizations  decreases. 

A  continuously  varying  medium  can  be  thought  of  as  a  series  of  thin  slabs 
of  different  n^.  Extraordinary  waves  are  produced  at  each  interface  and,  in  the 
limit,  the  coupling  region  may  be  thought  of  as  a  distributed  source  of  upgoing 
extraordinary  waves.  When  is  very  small,  the  extraordinary  waves  will  all 

have  nearly  the  same  phase  and  will  therefore  add  constructively, 

A  coupling  parameter,  xj/"  ,  is  defined  proportional  to  and  ^  in 

Po"Px 

accordance  with  the  above-mentioned  dependencies.  Using  the  fact  that 
PoPx  *  i  have; 


il-lnPo-^ 


p3  -  1  dz  2  dz  po  ♦  i 
From  the  Appleton-Hartree  equation,  the  polarization  p^  is  given  as; 

po  ‘  [  2Tu^"  yOT^xT^  } ' 


[20] 


[21] 
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Substituting  eqno  [21]  into  eqn«  [20]  we  obtain; 


r 


-iY^  Y 
4(U-X)^ 


dX 


[22] 


The  greater  the  coupling  parameter  ’'j/'  ,  the  larger  is  the  extraordinary  wave 
generated  in  the  coupling  region. 

At  X  ■  1,  the  denominator  of  eqn,  [22]  becomes  zero,  i.e.,  y  »  oo  (complete 
coupling)  at 


vi  2  4 

-4  Y7  +  Y^ 


(a) 


0; 


from  which  the  critical  collision  frequency  is  derived 


c  TcSsoT 


[23] 


Using  a  dipole  approximation  for  the  earth's  magnetic  field,  we  have; 


(1  +  3sin^0)^^^ 
r^ 


Using  the  above  relations. 


fn  -  2,8  Hp; 


tan  Of.  •  1/2  cot0, 
in  eqn,  [23]  we  have; 


which  is  used  in  Fig.  3  to  plot  vs,  the  geomagnetic  latitude  0, 
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Figure  3 


Transmission  coefficient  (T^j,) 

Let  us  define  €  such  that 

X  «  -1+62. 

JCrnf^ 

Close  to  the  ordinary  reflection  level,  X  1  (i,e»,  e  z  1)  and  we  can  utilize 
the  expansions 


^  *  ^0  *  ^  »  1  +  z 

from  which  it  follows  thats 


z 


The  transmission  coefficient  for  vertical  incidence  was  first  calculated  by 
Rydbeck  (1950) «  V,  L,  Ginzburg  (1961)  in  his  book,  pp,  515-559,  derives  it  in  three 


-19- 


different  ways.  We  will  avoid  repeating  the  derivations,  giving  only  the  final 
results. 


and 


where 


C 


_^3/2 


sin^a 
icosa  * 


Q  or  cosa  \/  y(Y  +  1) 

Po  *  o  - 5 — » 

^  2  (1  +  YcosO)^ 


c  »  3  X  10^®cm/sec. 


For  a  parabolic  layer,  where  N 


max 


[1  -( 


■Sir/ 


we  have 


f2  _  Ne^  _  *^max  ® 

*  Oc™  “  “31  m 


[1 


2 

1  » 
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An  ordinary  wave  of  frequency  f  will  be  reflected  at  X  »  1,  i,e.,  fj^  ■  f,  which  as 
seen  from  the  above  equation,  will  correspond  to  a  height 


For  a  parabolic  layer  therefore; 


f2 

Nmax 
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Thus,  for  V  ^  the  complete  expression  for  the  Z-mode  transmission  coefficient 
assuming  a  parabolic  layer  becomes: 


CrtAh  cos«VY(Y  ♦  1) 


“ITT 


(1  ♦  Y  cosO: )' 


f2  /£2 

*  ^^Nmax 
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When  ^oz  unity. 

Fig,  4  shows  ^  vs,  for  the  conditions  Ah  »  30  km,  0(  ■  13° 

C^geom,  lat,  65®),  fj^  ■  1,44  Me,  f  •  3,47  Me,  ^  ^c  “  0*23  Me 


Figure  4 
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From  FigSo  3  and  4  we  see  that  only  for  high  latitudes  (preferably  above 
70ON  or  below  70®S)  can  we  have  E-layer  coupling.  For  F-layer  coupling  where 
0.  even  higher  latitudes  are  required. 

For  example,  for  »  2f  »  4fj^,  Ah  a  40  km  and 

|^oz|^  “  ®*P  ("200  sin^)  which  means  that  at  cx  ci  9°,  i,e.,  geomagnetic  latitude 
0^  73°,  10"^,  At  Boston  where  0«  53°  and  cx  — 21°  |Tq2|^“  10"^^,  So  we 

see  that  we  can  have  vertical  coupling  only  at  high  geomagnetic  latitudes  in  the 
E-region  and  only  at  very  high  geomagnetic  latitudes  in  the  F-region, 

Oblique  incidence 

For  oblique  incidence  we  have  complete  transmission  when  sinB^  =  ^  ^  sino( 

(as  we  have  already  seen  from  eqn,  [19]),  with  a  collision  frequency  of  zero  as  for 

vertical  incidence,  when  OC  =  0.  If  collisions  are  present  they  will  broaden,  through 

coupling,  the  tolerance  in  the  angle  0^  allowing  a  cone  rather  than  a  single 

direction  of  propagation  in  the  Z-mode,  It  is  not  proper  to  calculate  this  cone, 

as  some  authors  have  done,  by  utilizing  the  critical  collision  frequency  for  vertical 
.  ,  ,  2XfHsin2o< 

incidence  «  2'co"o('""  *  formula  is  not  valid  for  oblique  incidence. 

The  problem  when  considering  oblique  incidence  is  not  the  penetration  of  the 
layer  by  the  Z-mode,  but  the  return  path.  One  has  to  assume  that  the  reflecting  layer 
is  rough  enough  to  produce  back  scattering,  so  that  some  extraordinary  wave  can  return 
along  the  incident  path  (parallel  to  the  magnetic  field) . 

The  observations  of  Ellis  support  this  theory  which  can  explain  the  existence 
of  the  Z-mode  at  moderate  latitudes,  where,  as  it  is  also  pointed  out  by  Budden  (1961), 
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the  explanation  through  vertical  coupling  would  be  very  difficult o  Some  doubts, 
however,  still  exist,  Ratcliffe  (1959), about  the  ability  of  the  reflection  layer 
to  produce  sufficient  back  scattering, 

Figo  5  demonstrates  the  generation  of  the  Z-mode  at  oblique  incidence  through 
back  scattering, 


It  is  quite  clear  that  the  larger  the  angle  0<  ,  i,e,,  the  lower  the  geomagnetic 
latitude,  the  more  difficult  it  is  to  have  sufficient  back  scattering  in  the  direction 
of  the  magnetic  field.  In  addition,  the  ionosonde  transmitters  used  were  primarily 
for  vertical  incidence,  so  that  the  intensity  of  the  obliquely  incident  wave  was  much 
weaker.  These  reasons  explain  why  the  Z-mode  has  not  been  observed  at  low 
latitudes  where  ©q  and  0<  must  be  quite  large. 
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Previous  work 

It  is  of  interest  to  give  an  account  of  the  work  done  to  date  on  the  Z-mode, 
discussing  at  the  same  time  some  of  the  ideas  put  forward  in  the  light  of  today's 
accepted  theories. 

To  Lo  Eckersley  (1933),  in  a  discussion  cf  the  ionosphere,  said  that  on  his 
wholly  visual  observation,  of  the  ionospheric  reflections  on  a  cathode  ray 
oscillograph,  he  noticed  that  "On  rare  occasions  there  are  very  definite  triplets, 
with  one  right  hand  and  two  left  hand  following,  all  sufficiently  close  to  ensure 
that  they  belong  to  the  same  system.  Triplets  have  been  observed  on  five  occasions 
during  the  last  five  months  on  wavelength  between  55  and  60  meters," 

G,  R,  Toshniwal  (1935)  reported  that  in  India,  working  on  a  single  frequency, 
they  observ  ed  mostly  after  sunset  triplets  lasting  for  not  more  than  a  minute, 

Toshniwal  still  uses  the  Lorentz  term  in  the  equations,  and  thinks  that  an 
evanescent  extraordinary  wave  partially  penetrates  and  reaches  the  other  reflection 
point,  which  Mary  Taylor  (1933)  had  suggested  as  a  possibility  if  the  electron 
density  varies  quite  rapidly  with  height, 

Leiv  Harang  (1936)  gave  the  first  records  of  ionograms  showing  a  Z-mode. 

This  is  also  called  triple  splitting  because  the  ionograms  have  the  form  shown  in 
Fig,  6, 
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Z  O  X 


-  f 


Figure  6 

He  observed  the  Z-mode  during  winter  at  Tromsb,  Norway,  only  during  the  day 
and  usually  around  noon  when  fQF2  was  high.  He  also  believes  he  is  observing  two 
extraordinary  and  one  ordinary  as  Mary  Taylor  had  suggested. 

Since  that  time  many  other  people  have  observed  the  Z-mode  and  have  published 
records  of  triple  splitting  from  different  parts  of  the  world. 

D.  H.  Meek  (1948)  reported  that  the  Z-mode  had  been  observed  in  Canada  since 
1943.  He  reports  that  it  was  often  seen  extending  from  the  E-region  through  the 
Fj  region  to  the  F2  region. 

From  a  sample  record  accompanying  his  report  he  reads; 

f2F2  ■  5,45  Me,  fQF2  «  6,10  Me,  -  6.95  Me. 

From  the  basic  equations  we  have  X  «  1  -  Y,  1,1+Y,  i.e,, 

fS  -  -  fxfH  ■ 
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solving  the  above  relations  for  we  obtain 


■  Cfx-fj)  -  2Cfx-fo)  -  — .  °  ■  2Cf„-fx)  ♦  -2^ 


which  yields  respectively 

■  1,50  Me,  1,60  Me,  1.38  Me, 

He  has  not  given  an  explanation  for  those  differences.  On  a  qualitative 
basis,  though,  they  can  be  explained  as  follows.  The  ordinary  and  extraordinary 
waves  are  deflected,  the  first  toward  the  North  Pole  where  the  critical  frequency 
is  lower,  and  the  second  toward  the  equator  where  the  critical  frequency  is  higher. 
Thus  (fjj-fg)  is  larger  than  it  would  have  been  if  both  were  reflected  at  the  same 
point,  and  thus  we  get  a  larger  value  of  fj^. 

The  critical  frequency  of  the  Z-mode,  however,  is  slightly  higher  because 
of  the  oblique  incidence,  i,e,,  (fz)obl,«^  ^^z^norm.  ThusCfo-^z) is  smaller  and 

COS  Oq 

gives  a  smaller  value  for  f^j.  Finally,  in  the  case  (fx”fz)»  both  are  increased  and 
thus  f|^  has  an  intermediate  value  probably  closer  to  its  actual  one. 

Triple  splitting  was  observed  also  by  Gordon  Newshead  (1948)  at  Hobart, 
Tasmania,  during  1947,  Only  39  clear  triple  splitting  records  out  of  5000  were 
identified,  occurring  mostly  between  the  hours  of  1700  and  2000  local  time, 

In  a  record  presented  he  gives: 


fjj  »  8,5  Me,  fp  -  7,6  Me,  f^  «  7  Me, 


He  states  that  is  higher  than  anticipated  and  claims  that  this  probably  can  be 
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accounted  for,  if  collisions  had  been  included  in  the  calculations.  Though 
collisions  might  play  some  role,  the  explanation  seems  basically  to  be  the  oblique 
incidence  of  the  Z-mode.  He  still  speculates  that  the  Z-mode  might  be  the  effect 
that  Mary  Taylor  suggested  in  1933, 

To  Lo  Eckersley  (1948),  commenting  on  Newshead's  observations,  refers  to  his 
original  observations  (1933)  where  two  of  the  three  received  waves  had  ordinary 
and  only  one  had  extraordinary  polarization.  If  Mary  Taylor's  scheme  was  the 
explanation,  we  would  have  two  waves  with  extraordinary  and  only  one  with  ordinary 
polarization.  He  concludes  that  the  explanation  of  the  Z-mode  is  the  vertical 
coupling  of  the  ordinary  and  the  extraordinary  waves  and  not  the  partial  penetration 
of  an  evanescent  extraordinary  wave  as  Mary  Taylor  had  suggested. 

The  Z-trace  polarization  was  accurately  measured  by  J.  E.  Hogarth  (1951)  at 
Fort  Chimo  in  Northern  Canada,  It  was  found  to  be  ordinary,  i.e.,  in  the  northern 
hemisphere,  anticlockwise  as  viewed  in  the  direction  of  propagation.  B.  Landmark 
(1952),  working  under  L,  Harang  in  TromsB,  Norway,  repeated  the  measurements  and 
his  results  agreed  completely  with  those  of  Hogarth.  Thus  it  was  clearly  established 
that  a  coupling  process  was  responsible  for  the  Z-mode. 

0.  E.  H.  Rydbeck  (1950)  derived  his  expression  for  the  transmission  coefficient 
due  to  coupling  at  vertical  incidence.  This  theory  can  explain  coupling  only  in  the 
E-region  and  only  for  high  latitudes  as  we  have  already  seen. 

Eckersley  (1950)  also  treated  the  coupling  phenomenon,  but  from  a  more  general 
point  of  view,  using  Rieman  surfaces. 
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Meanwhile,  James  C,  W,  Scott  (1950)  advocated  the  new  idea  that  the  Z-mode 

was  produced  by  propagation  at  oblique  incidence  along  the  magnetic  field  and  back 

scattering  from  a  rough  reflecting  layer.  He  implied  though  that  the  angle  of 

2  f  ° 

incidence  should  be  that  of  the  magnetic  field  and  that,  if  ^  ^  ^  ^  , 

we  have  longitudinal  propagation  which  holds  for  vertical  incidence  but  not 
necessarily  for  oblique  incidence,  as  is  the  case  here, 

G,  Ro  Ellis  (1953)  developed  more  carefully  the  theory  of  the  Z-mode  suggested 
by  Scott,  and  derived  the  correct  expression  for  the  angle  of  incidence,  6©,  for 
Z-mode  propagation;  ioC,, 


sinSQ  = 


0 


He  made  the  first  measurements  of  6^  in  Hobart,  Tasmania  (  Ot~18°)  using  two  spaced 
(120  fto)  loop  antennas  and  found  that  the  propagation  was  in  the  plane  of  the 
magnetic  meridian.  His  results  agreed  very  well  with  those  predicted  from  the  above 
formula. 

He  is  still  using  the  vertical  incidence  equations,  however,  and  therefore 
his  expressions  for  the  cone  angle  in  terms  of  the  critical  collision  frequency 
are  only  approximate. 

In  a  later  paper,  Ellis  (1956)  did  discuss  the  theory  of  the  Z-mode 
propagation  on  the  proper  basis  of  oblique  incidence  and  gave  measurements  of  the 
Z-mode  received  intensity  at  different  distances  from  the  transmitter.  The  results 
indicate  an  approximate  Gaussian  distribution. 
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In  the  same  article  he  pointed  out  that  the  hole  existing  in  the  ionosphere 
along  the  direction  of  the  magnetic  field  for  <^fo  could  be  of  use 

in  radio  astronomy.  For  a  difference  in  height  of  S  -km  we  have: 

fH-ffl  -  fH 


sine. 


Tpr-r 


sino( 


and 


sineQ 


i_slnc 


If'  ♦  f 


- 


TIT  .  f  ) 


therefore 


sin0Q-sin0o 


0O-9A  0„*e'  , 

2sin  cos  ■P^  o,  ^(eo-e;)coseo 


sS 


2R  fj^+f 


3  f 
TIT"  f 


in0. 


and  finally 


sino< 


^0  *  ^0-00  •  ♦  f  ‘a"®o« 

For  values  of  f  between  2fj^  and  3fj^  and  small  ©q  (^10°,  i,e,,  (X'«20°)  we  have 
@0—  1/6**  per  100  km.  If,  e.g,,  f^nax  *  fo  *  **  then  ^2  —  Operating 


at  t  3,5  Me  we  have 


For  a  parabolic  layer 


N  «  N 


max 


»-  max  -* 


^  y  I  therefore  ^  ■  1 1 


['  ■  (■=>'] 


or 


z 


7  ^max' 


That  is,  the  height  difference  of  the  two  coupling  regions,  above  and  below  the 
maximum  layer  is  about  half  the  thickness  of  the  parabolic  layer^  For  the  F-region, 
this  would  mean  more  than  100  km  which  would  imply  0^  0^2°o  Collisions  in  this 

region  are  very  few  making  uncertain  the  possibility  of  attaining  this  tolerance » 

Of  course  one  could  operate  much  closer  to  fp  to  reduce  6q,  but  attenuation  and 
irregularities  then  would  become  much  more  important, 

In  a  recent  article,  Ellis  (1962)  has  further  developed  his  previous  idea  for 
the  application  of  the  Z-mode  to  radio  astronomy  in  the  framework  of  present  day 
space  technology. 

He  pointed  out  that  a  satellite  orbiting  in  a  region  above  the  layer  of 

maximum  electron  density  could  receive  radio  waves  only  through  the  Z-mode  if 

it  operates  in  a  frequency  f  such  that  f ^  ^  f  ^  fjg  where  fj^  is  the  plasma  frequency 

of  the  environment  and  f,~  _JL_f  l/f2  -  filL  Receiving  radio  waves  only 

^  cos0o\lf'^ 

through  the  Z-mode  would  render  a  high  resolution  to  our  observations. 
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This  last  possibility  we  will  investigate  in  more  detail  in  the  second  part 
of  this  monographs 


Summary 

We  have  seen  the  basic  equations  for  propagation  of  electromagnetic  waves  in 

a  raagnetoionic  medium,  from  which  we  have  deduced  that  a  third  wave,  called  the 

Z-mode,  can  propagate  under  certain  conditions  beyond  the  height  at  which 
^  W  6  ^ 

f^  =  .  At  vertical  incidence,  the  Z-mode  will  be  produced  either  when  the 

%  m 

magnetic  field  is  almost  vertical  (high  geomagnetic  latitudes)  or  when  the  collision 

frequency  W  approaches  the  critical  collision  frequency  ioe,,  when 
.  2 

Vj,  ■  ZlTfjj  condition  is  met  only  in  the  lower  ionospheric  regions, 

provided  we  are  at  high  geomagnetic  latitudes,  so  that  is  not  too  large.  At 
mid-latitudes  this  coupling  process  is  extremely  weak  and  cannot  account  for  the 
phenomenon  of  the  Z-mode,  The  oblique  incidence  approach  must  therefore  be  adopted, 
For  this  case  we  have  seen  that  when  sinSo  »  ^  i  Z-mode  penetrates 

beyond  the  layer  X  ■  1  and  is  reflected  at  the  layer  X  »  (l  +  Y) ,  The  theory  assumes 
that  some  back  scattering  takes  place,  thus  allowing  part  of  the  wave  to  return 
through  the  same  path.  This  process  does  not  require  any  collisions  to  be  present. 

No  theory  for  oblique  incidence  has  yet  been  developed  to  give  the  transmission 
coefficient  for  angles  slightly  different  from  6^  »  sin’^^^y  sinO<^  or  for 
the  dependence  of  the  transmission  coefficient  on  the  collision  frequency.  One  can 


only  infer  from  the  existing  theory  for  vertical  incidence  that  the  cone  angle  about 
00  must  be  quite  small,  especially  for  the  very  low  collision  frequencies  that  prevail 
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in  the  upper  ionospheric  regions.  This  provides  the  Z-mode  with  strong  directivity 
and  renders  valuable  support  to  Ellis's  recent  proposal  for  its  use  in  obtaining 
high  resolution  satellite  radio  astronomical  observations  at  long  wavelengths. 
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